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From its inception, the Schwarzschild solution of Einstein’s equations has been the 
subject of controversy over the possibility of compressing matter within the gravitational 
radius 2GM. Given the accumulated evidence for supermassive compact objects ranging 
from a few IO^Mq to a few 10® Mq, the existence of black-hole-like objects is beyond 
doubt mu What does remain an issue is whether the Schwarzschild metric correctly 
describes the physics of the interior. Alternatives to classical black holes have been 
proposed with no singularities in the interior ismiainiiziEiiaiiniiiiiiiaiisi- 

The simplest model proposed for supermassive compact objects at the galactic 
centers is a self-gravitating degenerate fermion gas composed of, e.g., heavy sterile 
neutrinos pi m im IT31 iT^ However, this scenario cannot cover the whole mass range of 
supermassive black-hole candidates with a single sterile neutrino species na. 

Recently, Chapline et al uu put forth an interesting proposal based on analogies 
to condensed matter systems where the effective general relativity was an emergent 
phenomenon. Specihcally, they suggested that the sphere where the lapse function 
vanished marked a quantum phase transition, the lapse function increasing again at 
r < 2GM. As this required negative pressure, the authors of [Hj assumed the interior 
vacuum condensate to be described by de Sitter space with the equation of state p = —p. 

Subsequently, the idea of gravitational vacuum condensate, or ‘gravastar’, was taken 
up by Mazur and Mottola lai, replacing the horizon with a shell of stiff matter astride 
the surface at r = 2GM. Visser and Wiltshire HE] and recently Carter HE] also examined 
the stability of the gravastar using the Israel thin shell formalism [20]. Despite the fact 
that general relativity is an emergent phenomenon in string theory [21], the gravastar 
has met with a cool reception. Certainly, the assumption of a de Sitter interior presents 
a quandary: on the one hand, the quantum phase transition would suggest that the 
associated cosmological constant is a fundamental parameter; on the other hand, to 
accommodate the mass range of supermassive black hole candidates, it must vary over 
some six orders of magnitude. Thus, it seems prudent to explore other possibilities for 
the gravastar interior. 

One obvious extension of the de Sitter equation of state p + p = 0 would be 
an equation that satishes p + p < 0, thus violating the dominant energy condition. 
The fluid of which the equation of state violates the dominant energy condition has 
been dubbed phantom energy [221123 and has recently become a popular alternative to 
quintessence and to the cosmological constant [22] • The motivation for introducing the 
phantom energy is that the equation of state w =p/p < —1 produces a superaccelerated 
cosmological expansion which seems to be favored by the combined analysis of the 
CMB and super nova type I data [22] • However, a superaccelerated expansion may 
also be obtained without violating the dominant energy condition in scalar-tensor 
theories of gravity [22] and in models with variable gravitational constant [2Z] or variable 
cosmological constant [2H]- Some astrophysical aspects of the phantom have recently 
been discussed, such as phantom energy accretion [22] ; phantom energy wormholes n. 
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and a possible relation of the phantom tachyon model to supermassive black holes HQ. 

In this paper we consider a gravastar interior consisting of a self-gravitating scalar 
field described by a Born-Infeld type Lagrangian which yields the Chaplygin gas eqnation 
of state jSSES]- Hence, we look for static solutions of the self-gravitating Chaplygin gas. 
In particular, we consider static Chaplygin gas conhgurations in the phantom regime, 
i.e., when p-|-p < 0, and we show that these conhgurations could provide an alternative 
scenario for compact massive objects at galactic centers. 

The paper is organized as follows: In section |21 we introduce the basics of the 
model. In section El we investigate static solutions using Tolman-Oppenheimer-Volkoff 
equations. In section 0 we discuss a possible interpretation of galactic centers as Born- 
Infeld phantom gravastars. A stability analysis is given in section El and we conclude 
the paper by section El 


2. The model 


Consider the equation of state 

( 1 ) 

in the phantom regime, i.e., when 

P<VA. (2) 

Equation CD describes the Chaplygin gas which, for p > a/A, has attracted some 
attention as a dark energy candidate isaisni. Astrophysical objects made of the so- 
called generalized Chaplygin gas [311 have recently been discussed lH]. The generalized 
Chaplygin gas has also been exploited in the phantom regime |BS|. As we shall shortly 
demonstrate, static solutions to Einstein’s equations with matter described by CD with 
p)l cover the range of masses and radii required to fit the phenomenology of supermassive 
dark compact objects at the galactic centers. Moreover, Eq. CD yields the de Sitter 
gravastar solution in the limit when the central density of the static solution approaches 
the value \/~A. 

The Chaplygin gas equation of state CD with the condition p < \fA may be derived 
from the Dirac-Born-Infeld type Lagrangian 

-^^061 = —a/1 -|- X , (3) 

where 


X = . (4) 

Clearly, in the limit X —y 0 this Lagrangian becomes a free scalar held Lagrangian with 
a “wrong sign” kinetic term, hence a phantom||. A phantom Lagrangian of the type (E)) 

II It has recently been shown that quantum effects could lead to an effective dark-energy equation of 
state violating the dominant energy condition on cosmological scales even for a scalar field Lagrangian 
having the correct sign of the kinetic energy EZj. 
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has been proposed in the context of the superaccelerated expansion 
Lagrangian © leads to a perfect fluid with the four-velocity 

^ ,OL 

the pressure p = £dbi, and the density 


Note that the 


(5) 


P = 


< \fA. 


( 6 ) 


yiTx 

Clearly, the pressure and the density obey Eq. m and p < \/~A. It is worth noting that 
the Chaplygin gas cosmological model ia which the condition p > \fA holds, is 
described by the Lagrangian with a (—) sign in front of X. It is as though instead of 
the lapse function changing its sign, as it does in the Schwarzschild case for r < 2GM, 
the scalar kinetic energy changed its sign to become a phantom. 

Equivalent to £dbi of Eq. © is 

(T) 



as may be seen by eliminating 0^ through its equation of motion. Now recall that in 
four dimensions, a three-form held strength is dual to a scalar: 




,/3 


( 8 ) 


Then the content of the model is 


C 


1 

1202 


TT TT/lUa _ 1 



(9) 


Notably, this is also the content of low-energy string theory j2Ij in the Einstein frame 
where is the Kalb-Ramond held, 0^ corresponds to the string coupling and the 

dilaton kinetic term is neglected in lieu of a potential assumed to arise from non- 
perturbative ehects^. A potential similar to that in o and © has been considered in 
the context of polymer scaling and black holes ng. 


3. Static solutions 

Next, we proceed to solve Einstein’s equations for a static, spherically symmetric 
conhguration with an interior described by o and with a Schwarzschild exterior. Our 
approach is similar to that of Armendariz-Picon and Lim jlU]- However, in contrast 
to us, they consider the classes of Lagrangian with strictly space-like gradient of 0, 
i.e. with X < 0, so that in their case the quantity © cannot be interpreted as a four 
velocity. As a consequence, their energy-momentum tensor does not have a standard 
perfect huid form and the pressure is not isotropic. Here, we consider strictly time-like 
gradient of 0 with X > 0, but we allow phantom-like Lagrangians which violate the 
dominant energy condition. 

^ Strictly, the exterior should have = 6A and 0^ = q to be Schwarzschild. 





Born-Infeld Phantom Gravastars 


5 


For the static, spherically symmetric line element 

ds^ = dt^ - - - ' , -+ sin^ -d dip'^) 

1 - 2GM{r)/r ^ > 

with = diag (p, —p, —p, —p), Einstein’s equations become jH] 

= 47rr^ p, 

d - n c -^ + 4vr p 
^ ^ r(r-2GM) ’ 


( 10 ) 


( 11 ) 

( 12 ) 


while = 0 gives 

P' = -(p + p) j. (13) 

We focus on the equation of state m to close the system and we require the solution 
to (dD and m to be regular at r = 0. Since by rescaling t one may set ^(0) = 1, Eqs. 
m and (fT!^ yield 


= - 


a-pI 


Po ^ - p2 ■ 

Combining Eqs. o, m, and (0, one has 


P=G 1 


A 


'pM 


r(r 


2GM) J ■ 


(14) 


(15) 


In Figs. m m and 01 respectively, we exhibit the resulting p{r)/pQ, ^(r), and 
2GM.{r)/r for selected values of Pq/\/~A. The solutions depend essentially on the 
magnitude of po relative to \/~A. In the following we summarize the properties of three 
classes of solutions corresponding to whether po is larger, smaller, or equal to \/~A. 

i) For Po > -v/vd, the density p increases and the lapse function ^ decreases with r 
starting from the origin up to the black-hole horizon radius i?bh, where 2GAf(i?bh) = 
/?bh- In the limit po —oo, a limiting solution exists with a singular behavior 


p(r) ~ ( 


7A 


ISTrGr^' 


vl/3 


(16) 


near the origin. 

ii) For Po < \/~A, both p and ^ decrease with r up to the radius Rq where they 
vanish. At that point the pressure p blows up to — cxo owing to (HD- The enclosed mass 
M. is always less than r/(2G), never reaching the black-hole horizon, i.e., the radius 
where 2GM(r) = r. 

iii) For po = a/A, the density p remains constant equal to a/A up to the de Sitter 
radius i?ds = 2GM = 87rGR^gPo/3. Hence, the interior is de Sitter, precisely as in 
Chapline et ah jB]. The lapse function is given by 

1/2 


2 










Born-Infeld Phantom Gravastars 


6 



Figure 1. Density profile of the Chaplygin star for pojy/A = 1.2 (short dashed), 
0.98 (dotted), 0.9 (long dashed). The limiting singular solution with r~'^P behavior at 
small r is represented by the dot-dashed and the de Sitter gravastar by the solid line. 



Figure 2. Lapse function for various solutions as in Fig.^ 
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Figure 3. Enclosed mass divided by the radius of the star for various solutions as in 
Fig.m 


with 


-RdS = 


SttG 




( 18 ) 


As po from above or from below, solutions i) or ii), respectively, converge to 

iii) except at the endpoint. The lapse function in iii) joins the Schwarzschild solution 
outside 

^(0 = , (19) 

continuously, whereas in i) and ii) it happens discontinuously. 

As in the case of a de Sitter gravastar, in order to join our interior solutions to a 
Schwarzschild exterior at a spherical boundary of radius R, it is necessary to put a thin 
spherical shell'*' at the boundary with a surface density and a surface tension satisfying 
Israel’s junction conditions |2D|. In all the solutions discussed above the pressure is 
isotropic and does not vanish at the boundary. Hence, the pressure at the boundary 
must be compensated by a negative surface tension of the membrane. We postpone this 
issue for section El where we discuss the stability of the solution. 

It has recently been demonstrated that the joining can be made continuous without the presence of 
a thin shell for a gravastar made of the fluid with an anisotropic pressure m 
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4. Black holes at galactic centers 


Case ii), together with iii), is of particular interest as we would like to interpret the 
supermassive compact dark objects at the galactic centers in terms of phantom energy 
rather than in terms of a classical black hole. It is natural to assume that the most 
massive such object is described by the de Sitter gravastar, i.e., solution iii) (depicted 
by the solid line in Figs. 1, 2, and 3). If we identify the most massive black-hole 
candidate observed at the center of M87, with mass M^ax = 3 x with the 

de Sitter gravastar, then = O.TkeV"^, to be contrasted with the 10“^eV values 
wanted for cosmology jSB]- The radius of this object is i?ds equal to the Schwarzschild 
radius Clearly, solutions belonging to class ii), can fit all masses M < Mmax- 

However, for the phenomenology of supermassive galactic centers it is important to 
hnd, at least approximately, the mass-radius relationship for these solutions. This may 
be done in the low central density approximation, i.e., po ^ V^, which is similar to 
the Newtonian approximation but, in contrast to the Newtonian approximation, one 
cannot neglect the pressure term in Eq. m- Moreover, as may be easily shown, in 
this approximation Ai <C r^p, so that the pressure term becomes dominant. Next, 
neglecting 2GAt with respect to r, as in the usual Newtonian approximation, Eq. m 
simplihes to p' = AnGAr, with the solution 

..2 \ 



P = Po 

which gives a mass-radius relation 
M levr^ 


Rl = 


Po 


2nGA ’ 


Rl 


15 


GA = constant. 


( 20 ) 


( 21 ) 


The mass-radius relationship M (x. R^ which phantom gravastars obey, offers the 
prospect of unifying the description of all supermassive compact dark objects that have 
been observed at the galactic centers, as Born-Infeld phantom gravastars with masses 
ranging from Mmin = IO^Mq to M^ax = 3 x IO^Mq. Indeed, assuming that the most 
massive compact dark object, observed at the center of M87, is a Born-Infeld phantom 
gravastar near the black-hole limit, with i?max = 2G Mmax = 8.86 x lO^km = 8.21 Ihr, 
the compact dark object at the center of our Galaxy, with mass Mqq = 3 x IO^Mq, 
would have a radius Rqq = 2.06 Ihr if the scaling law m holds. This radius is well 
below the distances of the closest approach to Sgr A* which the stars SO-2 (i?min = 17 
Ihr = 123 AU, jlH]) and SO-16 (i?min = 8.32 Ihr = 60 AU (Hj) recently had and beyond 
which the Keplerian nature of the gravitational potential of Sgr A* is well established. 


5. Dynamical stability 

None of the solutions discussed in the preceding section will be stable unless there is 
a membrane, e.g. in the form of a thin shell, placed at the boundary with surface 
density a and surface tension 6 satisfying Israel’s junction conditions. Following Visser 
and Wiltshire ra we consider a dynamical thin shell allowed to move radially at the 
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boundary of the phantom gravastar and discuss under which conditions will a gravastar 
conhguration will be stable. A dynamical thin shell connecting two general static 
spherically symmetric spacetimes have also recently been considered SSI in a slightly 
different context. 

Israel’s junction conditions read 


/l'» _ 


= StiGS^ 


( 22 ) 


where = diag (a, 9, 6) is the surface stress energy, and [[/]] denotes the discontinuity 
in / across the boundary, i.e.. 


[[/(r)]] = lim {f{R + e) - f{R - e)) . 

e^O 

The tensor Kab is the extrinsic curvature dehned by 
Kab ~ 


(23) 

(24) 


where Uq is a spacelike unit vector orthogonal to the timelike boundary and hab is the 
induced metric on the shell. 

The angular components of the extrinsic curvature may be easily calculated from 
(121 yielding 


Kt = K^ = l{A + RY/\ (25) 

precisely as in UBI, where the dot denotes the derivative with respect to the proper time 
and where 


A 


2GM 


r 

The calculation of the time-time component 
easily done following Israel [20]. By making 
orthogonality from (121 it follows 


(26) 

Kl is slightly more involved and is most 
use of the Gauss normal coordinates and 


Kl 


u‘'u^Kab = 


(27) 


where is the four-velocity of a point on the shell. Its non-vanishing components for 
the metric m are 

1/2 




R 

aJ 


= R. 


Similarly, the non-vanishing components of Ua are 


nt = 


^R 


Hr = 


(A + i?2)l/2 


(28) 


(29) 


AV2 ’ A 

Using these expressions it may be shown that the four-acceleration u^u°‘-b satishes m 

1 -- -r ..b.M ( 30 ) 


TlaU U -^b 


(A + i?2)l/2 

A straightforward calculation yields 

1 


(i? + r>V). 


k: = 


(A + /?2)1/2 


f' 






V r / 



(31) 
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Then, from (P|) using m and m we find 


dvrcr 


1 

G 


\Ja + r^ 

R 


(32) 


47r(cr — 29) 


R^R + GM - AttGR^p AttrJa + R^ , 

- ;-- -^ + - (P + P) 

A + R^ ^ 


■ (33) 


To derive the stability condition a la Visser and Wiltshire ^H], it is now sufficient to 
replace IR by -2V{R) and R by -V'{R) in equations (1221) and (1221), where V{R) is a 
potential. The shell will be stable against small radial perturbation if there exists an 
equilibrium position R such that 

V{R) = 0; V\R) = 0; V"{R) > 0. (34) 

Then, by choosing a suitable potential V(R), equations (IH2|) and (I22|) dehne a parametric 
equation of state 9 = 9(a) for the shell. 

For the static shell and the metric m we hnd 


1 - 


R ) 


2GM(R)y^^ 


—AirGRa, 


and with the help of Eq. m we obtain 


M 1 - 


R ) 


= 47iR^(a -29), 


- M(i?) 


47ri?M\ / 2GM(R)\ 

I R ) 


(35) 


(36) 


where M is the total mass. Note that if the joining is affected at the point Rq where 
^ = 0 and p is inhnite (the point of naked singularity |l2]), the surface density a is 
hnite, whereas the surface tension 9 —> —oo. Hence, to avoid the singularity the shell 
must be placed at some R < Rq. 

For the Newtonian gravastars discussed in the preceding section, equations (1221) 
and o may be considerably simplihed. Using (1201), the approximation 2GM./R -C 1, 
and assuming that the boundary radius R is close to Rq, i.e.. 


Rq — R 

y = — < 1- 


R 


0 


we hnd 


a = 2 'kGAR, 


3 y 


VI - 2U ’ 


9 = - AnGARl^^= + t^GAR^^^ 


StiGRq y 


WI-2V 


’1-2U 


(37) 

(38) 

(39) 


In general, the potential U is a function of both Rq and y, depending on the chosen 
equation of state 9 = 9(a). The dynamical stability is achieved if for any Rq, V has a 
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minimum at some point y = 1 — R/R which should not depend strongly on Rq and at 
which V{y) = Q. 

Instead of postulating 9 = 9{a) we can choose a desirable potential and determine 
the equation of state a posteriori. Static stability (indifferent balance) is obtained by 
setting V = 0. Depending on the potential V{RQ,y), this may still be a complicated 
parametric equation of state. However, for the static shell, i.e. H = 0, the equation of 
state simplihes to 


9 = 


f ARq \ 


1/2 


[32nGa) 

This equation describes a 2-dimensional generalized anti-Chaplygin gas 
equation of state of the form p = —9 = A 'with a = 1/2. 


(40) 
with the 


6. Conclusions 

In conclusion, we have shown that replacing the de Sitter interior of the gravastar 
by a Born-Infeld phantom allows a wide range of gravastar mass and radii related to 
the central density, or equivalently to the velocity of the phantom scalar. We have 
demonstrated that if the constant H, as the only free parameter of the model, is hxed 
by assuming that the most massive galactic center object is the maximal (de Sitter) 
gravastar, then the model is able to explain all supermassive compact dark objects at 
the center of the galaxies. Furthermore, as demonstrated above, the phantom gravastar 
model can lay claim to a connection with low-energy string theory. 
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